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ABSTRACT 

Context. Classical Cepheids are excellent tracers of intermediate-mass stars, since their distances 
can be estimated with very high accuracy. In particular, they can be adopted to trace the chemical 
evolution of the Galactic disk. 

Aims. Homogeneous iron abundance measurements for 33 Galactic Cepheids located in the 
outer disk together with accurate distance determinations based on near-infrared photometry are 
adopted to constrain the Galactic iron gradient beyond 10 kpc. 

Methods. Iron abundances were determined using high resolution Cepheid spectra collected with 
three ditferent observational instruments: ESPaDOnS@CFHT, Narval@TBL and FEROS@2.2m 
ESO/MPG telescope. Cepheid distances were estimated using near-infrared (J, H, K-band) 
period-luminosity relations and data from SAAO and the 2MASS catalog. 
Results. The least squares solution over the entire data set indicates that the iron gradient in the 
Galactic disk presents a slope of -0.052±0.003 dexkpc"' in the 5-17 kpc range. However, the 
change of the iron abundance across the disk seems to be better described by a linear regime inside 
the solar circle and a flattening of the gradient toward the outer disk (beyond 10 kpc). In the latter 
region the iron gradient presents a shallower slope, i.e. -0.012±0.014 dex kpc"' . In the outer disk 
(10-12 kpc) we also found that Cepheids present an increase in the spread in iron abundance. 
Current evidence indicates that the spread in metallicity depends on the Galactocentric longi- 
tude. Finally, current data do not support the hypothesis of a discontinuity in the iron gradient at 
Galactocentric distances of 10-12 kpc. 

Conclusions. The occurrence of a spread in iron abundance as a function of the Galactocentric 
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longitude indicates that linear radial gradients should be cautiously treated to constrain the chem- 
ical evolution across the disk. 

Key words. Stars: abundances - Stars: supergiants - Galaxy: abundances - Galaxy: evolution 



1. Introduction 

Chemodynamical evolutionary models provide interesting predictions concerning the formation 
and evolution of the Milky Way. However, the plausibility and accuracy of these predictions need 
to be validated with observational data. The most commonly adopted observables are the star for- 
mation rate and the abundance gradients across the Galactic disk. 

Different tracers have been used to determine the Galactic gradients : HII regions, Open 
Clusters, (9/B-type stars and Cepheids. Compared to other tracers, Cepheids present several ad- 
vantages : i) they are well-known primary distance indicators; ii) they are bright enough to allow 
the study of the gradient over a large range of Galactocentric distances; Hi) they present a large set 
of well defined absorption lines, therefore, accurat e abundances of many elements can be provided. 

Abundance gradients were discovered first bv lSearld (1 197 lb in six external galaxies using HII 
regions. However, their occurrence in our Galaxy remained controversial as they are not easy to 
observe due to the fact w e are embedded in the G ala ctic disk. Ind eed, some studies were in favor 
of Galactic gradients, like lD'Odorico et aL ( 1976) or Janes ( 1 979 ) whereas others fou nd no corre- 



lation between distance and metalUcity (IClegg & BeU , 



1973 



Jennens & Heifer 



19751) . Today, the 



existence of Galactic abundance gradients seems widely accepted, even if some recent studies still 
report no evidence of such gradients, but the empirical determinations of their shapes and slopes 
are strongly debated. 

Considering a simple linear gradient the slopes for different elements can be explored with 
several stellar trac ers. The reader i n tereste d in exhaustive reviews con cerning Galactic gradi- 
ents is referred to lAndrievskv et al.l (l2002al) and IChiappini et al.l (1200 lb . By using HII regions, 
(Il996h found a slope of -0.02 dex kpc~', but other authors using the same trac- 



Vilchez & Esteban 



ers suggested slopes ran ging from -0.039 dex kpc ' (IDeharveng et al.L uOOOl) to -0.065 dex kpc ' 



dAfflerbach et al 



1997h. On the othe r hand, the use of B -type stars gives slopes ran ging from - 



0.042 dex kpc"' JPaflon et all |2004 to -0.07 dex kpc ' Joummersbach et al.lll998l) . The slopes 



based on planetary nebulae range from -0.05 dex kpc ' dCosta et al 



dMaciel et al 



1999b to the lack of a Galactic metallicity gradient dStanghellini et al 



2004 to -0.06 dex kpc 



2006). The 



slopes based on old Open Clusters still show a large sp read. By adopting a sample of 40 clusters 
distributed between the s olar circle and Rg - 14 kpc. 



Friel et al 



dex kpc . More recently. 



Carraro et al 



d2002l) found a slope of -0.06 
(12007b using new accurate metal abundances for f ive old 

J2OO2I) found 



Friel et al 



open clusters located in the outer disk together with the sample adopted by 
a much shallower global iron gradient, namely -0.018 dex kpc The global iron slopes based on 
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Cepheids are very homogeneous, dating back to the first estimates by iHarrisI ( 11981 



1984 . who 



found a sl ope of -0.07 dexkpc the more recent es t imates provide slopes ranging frorn -0.06 



dex kpc jAndrievskv et al 



2002c; 



to -0.07 dex kpc dLemasle et al 



Luck et al 



2003 



Andrievskv et al 



Paper I 



2004 



Luck et al 



2006) 



20071) . hereafterO 

As far as the shape of the gradient is concerned, current findings are even more controver- 
sial. In particular, the hypothesis of a linear gradient is still widely deba t ed. Se veral investiga- 



tions based on different stellar tracers — HII regions. 



Twarog et al Id 19971) 



Costa et al 



Carraro et al 



(l2007l) : Cepheids, 



Vilchez & Esteban 



Andrievskv et al. 



( 19961); open clusters. 



(120041) ; planetary nebulae. 



(120041) — indicate a flattening of the Galactic gradient b eyond 10-12 kpc. Th is flat- 



tening is well reproduced by chemodynamical evolutionary models dCescutti et al 



20071). More 



recently. I Yong et al.l d2006h brought forward a new feature: the flattening may occur with two dif- 
ferent basement values, the first one at -0.5 dex, a value more metal-poor than previous studies, and 
the second one at -0.8 dex. According to the quoted authors, the latter one might imply the possibil- 
ity of a merger event. On t he other hand, independ ent studies do not sho w evidence of a fla t tening 
toward the outer disk, like iRolleston et al.l d2000l) based on O/B stars or lDeharveng et alj d2000l) 
based on Hll regions. Mo reover, a discontinuit y in the metallicity gradient at Re ~ 10 - 12 kpc 
has also been suggested bv lTwarog et al.l dl997h according to open cluster photometry. Instead of a 
regular decrease from the inner to the outer disk, they indicate a 2 zone distribution with an abrupt 
discontinuity of — 0.2 dex at about 10-12 kpc and a s hallower gr a dient inside e ach zone. This hy 



pothes is was supported by 



Andrievskv et al 



Luck et al. 



d2003l) and by 



Andrievskv et al 



d2002cl) 

d2004l) . Such variety may be due to the fact that the different tracers adopted to estimate the gradient 
present only a handful of objects toward the outer disk, i.e. at Re > 12 kpc. 

In this paper, we investigate the shape of the gradient in the outer disk and focus our attention 
on the transition zone around 10-12 kpc. The paper is organized as follows: in §2 we describe the 
data and the reduction strategy adopted to estimate the iron abundances and the distances. In this 
section we also discuss the uncertainties affecting both abundances and distances. In §3 we discuss 
the new iron abundances and compare them with previous measurements. In §4, we present the 
new Galactic iron gradients and compare them with previous estimates based on different stellar 
tracers. §5 summarizes current findings. 



2. Observations and data reduction. 

2.1. Observations. 

High resolution spectroscopic obse rvations were perfo rmed in 2006 with ESPaDOnS at the CFHT 



(27 stars) and in 2007 with FEROS (iKaufer et al 



19991) at the 2.2m ESO/MPG telescope in La Silla 



(6 stars). Another target (BK Aur) already included in the ESPaDOnS sample was also observed 
in 2007 with Narval at the 2m Telescope Bernard Lyot (TBL) in Pic du Midi in Southwest France. 
ESPaDOnS has a resolving power of 81,000 in the spectroscopic "star-only" mode (only the light 
from the star goes through the instrument) and an has accessible wavelength range of 370-1050 
nm. FEROS has a resolving power of 48,000 and an accessible wavelength range of 370-920 nm, 
while Narval is a copy of ESPaDOnS adapted to the specifics of the TBL, whose resolving power 
is 75,000 in the spectroscopic "star-only" mode and an accessible wavelength range of 370-1050 
nm. The spectra were reduced either using the FEROS package within MIDAS or using the Libre- 
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ESpRIT software dDonati et al. 



1997 



20061). Relevant information concerning the observations 



and the signal-to-noise (S/N) ratio at 600 nm are listed in Table [T] Phases were cal culated with 
period s and epochs from the GCVS (General Catalogue of Variable Stars) described in lSamus et al 



( 120041) . except for 8 stars to take into account changes in the period. In the case of AO Aur, RZ 
Gem, SV Mon, UY Mon, RS Ori, GQ Ori, the y are fromberdnikov & Ignatoval J2OO0I) while in 
the case of AD Gem and CV Mon they are from 



Szabadosl(ll99ll) . 



Table 1. Observations: date, phase, spectrograph, exposure time and S/N ratio at 600 nm. 



Object 


date 


Julian date 
d 


phase 


spectrograph 


exp. time 

s 


S/N 
(600nm) 


AO Aur 


14/02/06 


2453781.910 


0,632 


ESPaDOnS 


1500 


89 


AX Aur 


16/02/06 


2453783.775 


0.619 


ESPaDOnS 


2*2700 


113 


BK Aur 


14/02/06 


2453781.974 


0.150 


ESPaDOnS 


700 


101 




31/10/07 


2454405.736 


0.096 


NARVAL 


1800 


101 


SYAur 


15/02/06 


2453782.928 


0.809 


ESPaDOnS 


500 


125 


Y Aur 


16/02/06 


2453783.911 


0.983 


ESPaDOnS 


800 


142 


YZ Aur 


15/02/06 


2453782.712 


0.789 


ESPaDOnS 


1600 


145 


AOCMa 


31/03/07 


2454191.099 


0.241 


FEROS 


2700 


107 


AA Gem 


15/02/06 


2453782.818 


0.749 


ESPaDOnS 


900 


128 


AD Gem 


15/02/06 


2453782.971 


0.131 


ESPaDOnS 


900 


126 


RZGem 


15/02/06 


2453782.804 


0.768 


ESPaDOnS 


800 


91 


BE Mon 


15/02/06 


2453782.886 


0.409 


ESPaDOnS 


1600 


120 


BVMon 


16/02/06 


2453783.950 


0.851 


ESPaDOnS 


1800 


79 


CV Mon 


16/02/06 


2453783.928 


0.084 


ESPaDOnS 


1500 


152 


EKMon 


15/02/06 


2453782.838 


0.809 


ESPaDOnS 


1800 


103 


SV Mon 


15/02/06 


2453782.955 


0.663 


ESPaDOnS 


300 


119 


TW Mon 


31/03/07 


2454191.074 


0.664 


FEROS 


2*1800 


80 


TX Mon 


14/02/06 


2453781.932 


0.998 


ESPaDOnS 


1700 


111 


TYMon 


15/02/06 


2453782.737 


0.537 


ESPaDOnS 


1800 


72 


TZMon 


14/02/06 


2453781.889 


0.896 


ESPaDOnS 


1400 


106 


UY Mon 


16/02/06 


2453783.967 


0.065 


ESPaDOnS 


600 


107 


V495 Mon 


15/02/06 


2453782.911 


0.074 


ESPaDOnS 


1800 


46 


V508 Mon 


14/02/06 


2453781.954 


0.101 


ESPaDOnS 


1700 


130 


VSlOMon 


31/03/07 


2454191.029 


0.975 


FEROS 


2*1800 


98 


WW Mon 


16/02/06 


2453783.844 


0.884 


ESPaDOnS 


2*2700 


104 


XX Mon 


15/02/06 


2453782.863 


0.584 


ESPaDOnS 


1800 


58 


CS Ori 


15/02/06 


2453782.761 


0.637 


ESPaDOnS 


1800 


78 


GQOri 


15/02/06 


2453782.961 


0.836 


ESPaDOnS 


300 


125 


RS Ori 


15/02/06 


2453782.981 


0.709 


ESPaDOnS 


350 


118 


HWPup 


29/03/07 


2454189.081 


0.906 


FEROS 


4*1800 


93 


AVTau 


16/02/06 


2453783.722 


0.701 


ESPaDOnS 


2*1800 


73 


STTau 


15/02/06 


2453782.949 


0.696 


ESPaDOnS 


300 


120 


DD Vel 


26/03/07 


2454186.590 


0.329 


FEROS 


4*1800 


127 


EZ Vel 


26/03/07 


2454186.495 


0.106 


FEROS 


4*1800 


155 




30/03/07 


2454190.204 


0.214 


FEROS 


2*1800 


129 
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2.2. Line list 



5 



We started from the iron line list provided by lRomaniello et alJ (l2008h and used in lPaper J . This list 
was extended by following the method devised by the quoted authors to match the wider spectral 
range of ESPaDOnS and Narval. The atomic properties (oscillator strength, excitation potential) 
Usted in VALD (Vienna Atomic Lines Database) were adopted for all the selected hnes. 



2.3. Equivalent widths 



Measurements of equivalent widths ( EW) are mad e with a semi-interactive code based on genetic 
algorithms from lCharbonneaul d 1 995b (see lPaperj for details). The analysis takes into account only 
weak, non blended lines so that all of them can be fitted by a Gaussian. While strong lines are 
systematically discarded, weak and asymmetric lines are removed after visual i nspection. No pe- 



Andrievskv et al 



culiar treatment was applied to weakly asymmetric lines, like the one proposed by 
( l2005b . To ensure the reliability of our EWs, we compared them to EWs obtained by direct integra- 
tion of the lines in the case of RZ Gem. As no systematics were found between direct integration 
and Gaussian fit, we assumed the asymmetry was weak enough so that the lines could still be fitted 
by Gaussians. The number of lines finally used is on average 80-130 lines for Fel and 15-25 lines 
for Fell. 



2.4. Effective temperature estimates 



The determination of an accurate temperatur e is a k ey point in the abundance determina 
tion. Following the same approach adopted in iPaper j. the stellar effective temperatures, T^ff 



Kovtvukh & Gorlova 



(120001). 



were estimated using the method of line depth ratios described in 
This approach has the ad ditional advantage of being independent of interstellar reddening. 
Kovtvukh & Gorloval (120001) proposed 32 analytical relations to derive Tj,// from line depth ratios. 
After measuring the line depths and calculating associated ratios, we obtaine d 32 temperatures 
whose mean value gave us the stellar temperatures listed in Table |2] Recently iKovtvukhl (l2007h 
increased to 131 the set of relations for the determination of supergiant temperature. We double- 
checked our temperatures using these new relations (Kovtyukh, private communication) and we 
found a very good agreement between the two independent estimates. 



2.5. Atmospheric parameters 

The surface gravity, logg, and the micro-turbulent velocity, v,, were determined by imposing an 
ionization balance between Fel and Fell with the help of the curve of growth. Lines of the same 
element in different ionization states should give the same abundance value. Iterations on logg and 
V, are, therefore, repeated until Fel and Fell adjust to the same curve of growth. This tool also 
provides an independent check for the Te/f value, since both high and low Xe.^ values properly fit 
the curve of growth. Once abundances were measured, atmospheric parameters are validated by 
checking that the Fel abundances depend neither on line strength nor on the excitation potential. 
The atmospheric parameters for the Cepheids in our sample are listed in Table |2] 
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Table 2. Atmospheric parameters adopted to compute iron abundances. 



Object 


phase 




logs 


V, 


[Fe/H] 






K 


dex 


km s"' 


dex 


AO Aui- 


0,632 


5450 


0.8 


3.0 


-0.41 


AX Am- 


0.619 


5635 


0.8 


2.7 


-0.22 


BK Aui- 


0.150 


5970 


1.3 


3.0 


-0.10 




0.096 


6160 


1.4 


3.4 


-0.04 


SY Aur 


0.809 


6220 


1.1 


3.0 


-0.07 


Y Aur 


0.983 


6300 


1.5 


3.0 


-0.26 


YZ Aur 


0.789 


5715 


1.0 


3.5 


-0.33 


AOCMa 


0.241 


6065 


1.7 


3.7 


-0.04 


AA Gem 


0.749 


5500 


1.0 


4.7 


-0.35 


AD Gem 


0.131 


6130 


1.6 


3.0 


-0.19 


RZ Gem 


0.768 


5500 


0.8 


3.3 


-0.44 


BEMon 


0.409 


5865 


1.4 


3.0 


-0.07 


BVMon 


0.851 


6080 


1.9 


3.2 


-0.10 


CV Mon 


0.084 


6350 


2.0 


3.1 


-0.10 


EK Mon 


0.809 


5700 


1.3 


3.4 


-0.05 


SV Mon 


0.663 


4900 


0.5 


3.4 


-0.10 


TW Mon 


0.664 


5640 


1.3 


3.4 


-0.15 


TXMon 


0.998 


5945 


1.0 


3.0 


-0.12 


TYMon 


0.537 


5555 


1.2 


3.0 


-0.15 


TZ Mon 


0.896 


6020 


1.3 


3.0 


-0.04 


UY Mon 


0.065 


6100 


1.5 


3.0 


-0.33 


V495 Mon 


0.074 


6090 


1.6 


3.2 


-0.17 


V508 Mon 


0.101 


6100 


1.5 


3.2 


-0.25 


VSlOMon 


0.975 


5390 


0.8 


3.3 


-0.12 


WW Mon 


0.884 


6650 


1.4 


3.0 


-0.32 


XX Mon 


0.584 


5520 


0.9 


3.0 


-0.18 


CSOri 


0.637 


6350 


2.0 


3.1 


-0.19 


GQ Ori 


0.836 


4960 


2.6 


0.8 


0.11 


RS Ori 


0.709 


5475 


1.3 


3.0 


-0.14 


HWPup 


0.906 


5300 


0.6 


3.1 


-0.28 


AV Tau 


0.701 


5750 


1.3 


3.0 


-0.17 


STTau 


0.696 


5700 


1.4 


3.2 


-0.14 


DD Vel 


0.329 


5600 


1.1 


3.2 


-0.35 


EZ Vel 


0.106 


6380 


1.5 


3.6 


-0.01 




0.214 


5455 


0.6 


3.6 


-0.01 



2.6. Abundance determinations 



Atmosphe ric parameters are t hen us ed as input for stellar atmosphere models. We used the MARCS 



models of 



Edvardsson et al 



( Il993h which are based on the following assumptions: plane-parallel 
stratification, hydrostatic equilibrium and LTE. The abundance determination codes adjust abun- 
dances until a good match between predicted and observed equivalent widths is obtained. The final 
abundance of a star is estimated as the mean value of the abundances determined for each fine. The 
high S/N ratio of the spectra and the use of a substantial number of clean lines provide the oppor- 
tunity to reduce the intrinsic errors on iron abundance t o ~ 0.12 dex. Note that in this investigation, 
we adopted the solar chemical abundances provided bylj 



Grevesse et al. 



(119961) . 
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7 



The very first source of errors in abundance determination comes from the extraction of data from 
the spectra. The equivalent width (EW) and the continuum placement must be accurately measured. 
Thus, we limited our analysis to symmetric and unblended lines, avoiding the most crowded parts 
of the spectra. We retained only relatively weak lines (EW < 120 mA) for the abundance calculation 
and checked that they could be fitted by a Gaussian. 

The crucial point in the abundance determination is to obtain an accurate temperature as it 
strongly aff'ects the line stren gth. As mentioned above, we used the method of line depth ratios from 



Kovtvukh & Gorloval (l2000h . according to which the uncertainties on absolute temperatures are at 



most X 80 K. Moreover, the curves of growth provided an independent check of the temperature 
determination as lines with high and low Xex values must both fit well the curve of growth. This 
ensures us that our intrinsic dispersion in temperature is not more than 100 K. 

In order to test how the uncertainties aff'ect the final abundance result, we computed abundances 
by adopting over or under-estimated values of the atmospheric parameters. As expected, only the 
errors on temperature have a noticeable eff'ect: over or under-estimating the efi'ective temperature of 
100 K implies a diff'erence of about + 0. 1 dex in the abundance determination. EiTors on the surface 
gravity of Alogg = +0.3 dex and on the micro-turbulent velocity of Av, = +0.5 km/sec have modest 
effects as the differences are only ± 0.03 and + 0.05 dex respectively. The sum in quadrature of 
these uncertainties on the atmospheric parameters gives uncertainties on the abundances of ~ 0.12 
dex. 

Finally, the LTE assumption might not be appropriate for supergiant, variable stars. The adopted 
stellar atmosphere models could then lead to s ystematic error s . How ever, the use of relatively weak 
lines should minimize this effect. Moreover, iFrv & Carnevi (Il997h have shown that a canonical 



spectroscopic analysis, using LTE atmosphere models, gives reliable abundances for Cepheids 
dwarfs and Ceph eids located in th e same open clusters are found to have the same metallicity, 
within the errors. 



Yong et al 



(120061) also tested the influence of NLTE effects and showed that the 



surface gravities derived from a classical approach are robust. 



2.8. Distances 



Individual Cepheid distances were estimated using near-infrared (J, H, ^T-band) photometry either 
collected from SAAO or from the 2MASS catalog. The mean magnitude of the Cepheids based 



Soszvnski et al 



on 2M ASS photometry was estimated using the template light curves provided by 
(l2005b togeth er with the V-band amplitude and the epoch of maximum available in the literature. 
The data from 



data from 



Barnes et al. 



Schechter et al. 



(119971) were also used for ST Tau and AD Gem while for HW Pup, the 



(119921) were used together with 2MASS data. In the case of UY Mon and 



AO CMa that are first overtone pulsators, it is not possible to use the empirical template, since it is 
only available for fundamental mode pulsators. The mean near-infrared magnitude of these objects 
is only based on the 2MASS measurement. Their period was fundamentalized as explained below. 
Note that the SAAO magnit udes were transfor med into the 2MASS photometric system using the 
transformations provided by[ 



Koenetal 



(12007D 



We adopted the individual Cepheid reddening s provided either by 



Lanev & Caldwell 



listed in the Fernie database dFernie et al 



(l2007l) or 



19951) . These two sets of reddenings are not homoge- 
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Fernie et al 



neous, as recent results p rovided by Lanev & Caldwelll (l2007b are corrected for metallicity effects 
whereas reddenings from 



([19951) do not take into account metallicity effects. The rel- 
ative absorption in the V-ba nd was estimated us ing Ay = 3.1 E(B-V) with Ay = 0.28 Ay, Ah - 



0.19 Av, and A,^ = 0.11 Ay dCardelh et al 



1989h . 



The Cep heid distances were e stimated using the near-infrared Period-Luminosity (PL) relations 
provided bv lPersson et al.l (120041) and by ass uming an LMC distanc e modulus of 18.50, a conser- 



vative value adopted by the HST Key Project ( Freedman et al. 
distance modulus, please see 



Romaniello et al 



200 ih . For a discussion on the LMC 



(120081) and references therein. Note that the quoted 
PL relations were also transfo rmed into the 2MASS photometric system using the transformations 
provided by ICarpenten (1200 ih . In order to estimate the distance of first overtone Cepheids their 



periods were fundamentalized, i.e. we added 0.127 to their logarithmic period. 



Using these PL relations from 



Persson et al. 



(120041) for Galactic Cepheids is justified although 



they are based on a large and homogeneous sample of LMC Cepheids. CuiTent empirical and theo- 



dent on metallicity dBono et al. 



1999 


Persson et al. 


2004 


Fouque et al. 


2007 


Romaniello et al. 



20081). Moreover, the near-infrared PL relations present an intrinsic dispersion that is significantly 
smaller than the dispersion of the optical PL relations. The difference is due to the fact that the 
former ones are, at fixed period, marginally affected by the intrinsic width in temperature of the 
instability strip. Finally, Cepheid distances based on near-infrared PL relations are more accurate 
than optical ones because they are much less affected by uncertainties in reddening estimates. 

The heliocentric distance of the stars adopted in the following is the mean value between 
the three distances in /, H, and K bands to limit random errors. The Galactocentric distance is 
calculated with the classic al formula assuming a Galactocent ric dista nce of 8.5 kpc for the sun 



dFeast & Whitelockl 



19971) . in order to enable comparisons with lPaper II This determination (whose 



error bars are +0.5 kpc) was based on the analysis of the Hipparcos proper motio n of 220 Cepheids . 



Recen t studies give Gal actocentric distances f or the sun between 7.62+0.32 kpc dEisenhauer et al 



2005b and 8.8+0.3 kpc dCoUinge et al 



Center, see 



Groenewegen et al 



d2008l) 



20061) . For a recent review of the distance to the Galactic 



The distances based on near-infrared photometry are given in the second column of Table [3] 
Galactocentric distances based on V-band photometry for which we adopted the heli ocentric dis 
tances given in the Fernie database are given only as a guideline in the first column (see 



Fernie et al 



d 1995b for details). Indeed these distances were estimated using a V-band period-luminosity re- 
lation but this relation is not corrected for metallicity effects and was derived using individual 
reddenings determined from a period-colour relation which doesn't acco unt for the metallicity . 



whereas these per iod-colour relations are likely to be metallicity-dependent dLanev & Stobie , 



1994 



Bono et al 



1999) 



As expected, the two different sets of distance determinations agree quite well near the solar circle, 
where the metallicities are in general near to the solar value. However, for lower metallicity objects 
located in the outer disk, the discrepancy increases and reaches respectively 2.8%, 3.5%, 3.4% and 
4.8% in the case of V510 Mon, TW Mon, HW Pup and YZ Aur, 4 stars among the 5 more distant 
Cepheids in our sample. 
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Table 3. Cepheid distances: 

1st column: Galactocentric distance from V band photometry. 
2nd column: Galactocentric distance from IR photometry. 
3rd coluimi: Hehocentric distance from IR photometry. 



Object 


Galactocentric 


Galactocentric 


Heliocentric 




distance (V) 


distance (IR) 


distance (IR) 




kpc 


kpc 


kpc 


AO Aur 


12.60 


12.43 


3.940 


AX Aur 


12.79 


12.75 


4.256 


BK Aur 


10.72 


10.86 


2.500 


SY Aur 


10.71 


10.92 


2.486 


Y Aur 


10.27 


10.42 


1.969 


YZ Aur 


13.30 


12.71 


4.278 


AO CMa 


1 1.32 


11.10 


4.068 


AA Gem 


12.39 


12.31 


3.820 


AD Gem 


11.20 


11.37 


2.939 


RZ Gem 


10.57 


10.66 


2.158 


BE Mon 


10.02 


10.11 


1.742 


BV Mon 


10.89 


11.19 


2.986 


CV Mon 


10.22 


10.09 


1.774 


EK Mon 


10.84 


10.71 


2.575 


SV Mon 


10.94 


10.82 


2.514 


TW Mon 


14,47 


14.00 


6.079 


TX Mon 


12.56 


12.34 


4.344 


TY Mon 


11.90 


11.94 


3.890 


TZ Mon 


12.22 


12.01 


3.982 


UYMon 


10.74 


10.37 


2.038 


V495 Mon 


12.79 


12.65 


4.679 


V508 Mon 


11.38 


11.43 


3.222 


V510 Mon 


13.73 


13.38 


5.337 


WW Mon 


13.68 


13.52 


5.274 


XX Mon 


12.00 


12.33 


/I ^ 1 T 

4.01 / 


CS Ori 


12.42 


12.56 


4.214 


GQOri 


10.91 


10.83 


2.469 


RSOri 


10.00 


10.02 


1.575 


HWPup 


13.28 


13.77 


7.795 


AVTau 


11.51 


11.54 


3.043 


STTau 


9.50 


9.52 


1.044 


DD Vel 




10.69 


6.710 


EZ Vel 




12.20 


9.522 



3. Results : Iron abundances 

A large fraction of Cepheids in our sample present slightly sub-solar iron abundances, i. e. the 
metal content expected for Cepheids located in the outer disk. Current results are in very good 
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agreement with previous stud ies, mainly t he ser i es of papers from Andrievsky a nd collaborators 



( Andrievsky et al. 




2OO2J 


I 


2005 Jb; 


Luck et al 




2006). 1 



Luck et al 



2003 



Andrievsky et al 



2004 



2005 



Kovtvukh et al 



dances. On the other hand, for six stars the differences are ~ 0. 1 dex and only for four stars (UY 
Mon, BK Aur, RZ Gem and AO Aur) is the difference higher and approaches 0.3 dex. 

The differences could arise from the methods employed to determine abundances, as 
Andrievsky and collaborators used a modified version of the standard spectroscopic analy- 
sis in order to avoid a possible dependence of Fel lines on NLTE effects, as explained in 



Kovtvukh & Andrievskvl ( 119991) . In their method, surface gravities and microturbulent velocities 
are derived from Fell hues instead of Fel lines. However, the differences are not correlated with 
specific phases, atmospheric para meters or metalhcities. Another di fference between the studies is 



the value of atomic parameters, as iKovtvukh & Andrievskvl (119991) derived their logg through an 



inverted analysis of the solar spectrum. However, logg for Fel are very similar for lines in com- 
mon. Our list of Fell lines is about twice as large due to the wider spectral range, and there are 
only mild differencies for lines in common, except for 2 lines where they reach x 0.2 dex. Thus 
differences in atomic parameters cannot be the origin of the discrepancies in iron abundances. 

The curves of growth for Fel (Fig.[T]i and for Fell (Fig.|2]l are presented below in the case of RZ 
Gem. The good agreement between predicted curves of growth and measurements for both Fel and 
Fell as well as the test performed on EWs measurements (see §2.3) indicate that our abundances 
measurements are robust. The same outcome applies for the other discrepant Cepheids. However, 
these stars deserve further investigation as they are all located in the outer disk, which is less well 
sampled than the solar neighborhood. 




-7 -6 -5 -4 -3 

log (ttj, * a! * r) 

Fig. 1. Observed curve of growth for Fel in RZ Gem. The full line is the theoretical curve of 
growth for a typical line (A-5QQQA, Xex=^)- The atmospheric parameters adopted for this star are 
Teff-55QQ K, log g=0.8, v,=3.3 km/s, and [Fe/H]=-0.44 dex. Squares represent lines with Xex < 
1.5 crosses lines with 1.5 <Xex < 3.0 and pluses hues withxex > 3.0. 
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log t « D 



Fig. 2. Same as Fig.[T]but for Fell. 



A few stars have been investigated in many studies and the results are quite often homogeneous 
(see, e. g., SV Mon). Only three stars of our sample ( YZ A ur, WW Mon and HW Pup) are in 



the values provided bv lYong et al 



Luck et al 



common with both Andrievsky's series and Yong et al. 
very similar to the analysis performed b y 
well with those publis hed in 



Yon 



g et al 



( 20061) . Although our classical approach is 



( 20061) . current iron abundances agree quite 



(l2006l) . 



(120061) ■ We found metal abundances that are higher than 



This difference cannot be explained by the different approaches adopted in the two investi- 
gations to estimate the effective temperature. Indeed, independent support of current Tf// values 
is provided by the evidence that both high and low Xex lines properly fit the curve of growth. 
Moreover, we found that the two investigations adopted different ^/-values for the Fell lines, but 
the difference is on average smaller than 0. 1 dex and the ensuing difference in the estimated sur- 
face gravities is negligible. Therefore, this cannot account for the differences in iron abundances. 
Although our spectra have higher spectral resolution (48000 -81000 versus 280 00) and higher S/N 
ratio (145, 104, 93 versus 49, 60, 49) than the spectra fror nYong et al.l (120061) . the discrepancies 
correlate neither with the S/N ratio nor with the atmospheric parameters. We cannot reach any firm 
conclusion regarding this discrepancy, since there are only three objects in common in the two 
Cepheid samples. 



4. Results : Galactic abundance gradients 

4.1. The iron gradient along tlie wliole Galactic disk 

The current sample inc ludes 33 Cepheids an d we added to this new data the 30 Cepheids of our 
previous investigation dLemasle et al.l l2007h . The 63 stars of this sample are located between 8 
and 15 kpc from the Galactic center (see Fig. [3]l. We investigated at first the working hypothesis 
of a linear gradient and calculated the iron Galactic gradient using the entire Cepheid sample. We 
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Andrievskv et al 



Luck et al 



Luck et al 



(2003), (4) 



(2006), (6) 



Kovtvukh et al 



Ponll(ll997b . 



(l2002al). (2) I Andrievskv et al 



Andrievskv et al. 



Yong et al 



(l2005ah . (8) 



oooel). 



(|2004 



(l2002cl) . 



Ta ble 4. Cepheid iro n abund ance and comparisons with pre vious results: 

(1) 
(3) 
(5) 
(7) 
(9) 



Kovtvukh et al 



(l2005bl) . 



Object 


[Fe/H] 


[Fe/H] : Previous results 




dex 


dex 


AO Aur 


-0.41 


-0.14^ 


AX Aur 


-0.22 




BKAur 


-0.10 


0.17^ 




-0.04 




SY Aur 


-0.07 


-0.02** 


Y Aur 


-0.26 


-0.23' 


YZ Aur 


-0.33 


-0.07', -0.37^ -0.64^ -036^ 


AOCMa 


-0.04 


-0.02' 


AA Gem 


-0.35 


-0.24^* 


AD Gem 


-0.19 


-0.19^ 


RZ Gem 


-0.44 


-0.12'' 


BE Mon 


-0.07 




BV Mon 


-0.10 




CV Mon 


-0.10 


-0.03' 


EK Mon 


-0.05 


-O.P 


SV Mon 


-0.10 


-0.03', 0.00^-0.04^ -0.02^ 


TWMon 


-0.15 


-0.24^ -0.20" 


TXMon 


-0.12 


-0.142 


TYMon 


-0.15 




TZ Mon 


-0.04 


-0.032, -0.12^ -0.24' 


UY Mon 


-0.33 


-O.OS'' 


V495 Mon 


-0.17 


-0.262,-0.16' 


V508 Mon 


-0.25 


-0.252 


V510Mon 


-0.12 


-0.1 9^ -0.20'' 


WW Mon 


-0.32 


-0.29\ -0.25', -0.59* 


XX Mon 


-0.18 


-0.182, -0.08^ 


CS Ori 


-0.19 


-0.262 


GQ Ori 


0.11 


-0.03', 0.06' 


RS Ori 


-0.14 


-0.l',-0.023 


HWPup 


-0.28 


-0.292, -0.2^, -0.40'' 


AVTau 


-0.17 




STTau 


-0.14 


-0.05' 


DD Vel 


-0.35 




EZ Vel 


-0.01 





found a slope of -0.023+0.007 dexkpc ' that is shallower than previous studies. However, this 
low value can be easily explained by the lack of targets inside the Solar circle : in this region, the 
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metallicities are higher an d indeed, inner Cepheids approach abundances of the order of +0.3 dex 



( Andrievskv et al. 



2002bl) . The inclusion of these Cepheids will then dramatically change the slope 



of the gradient. Moreover and even more importantly, data plotted in Fig.[3]do not show evidence 
of a gap over the Galactocentric distances covered by the current Cepheid sample. 



SLOPE: -0.025 dex/kpc 



0.5 



0.0 



-0.5 



-1 .0 



a 10 12 

Galactocentric distance ( kpc ) 



14 



Fig. 3. Galactic iron abundance gradient. The solid line shows the linear regression over the entire 
sample. Characteristic error bars are plotted in the lower right corner 



In order to improve the sampling across the Galact ic disk, we added to the se 63 Cepheids 



Romaniello et al. 



(l2008h . For the entire 



115 stars from Andrievsky's sample and 10 stars from 
Cepheid sample, we have both metallicities and infrared photometry. Among them, 38 Cepheids 
have Galactocentric distances between 10 and 12 kpc and 27 are even more distant. Individual 
distances and iron abundances are listed in Table 14.11 A map (as seen from above the Galaxy) 
showing the spatial distribution of the Cepheids is displayed in Fig.|4] The x-axis contains both the 
center of the Galaxy at (0,0) and the Sun at (8.5,0). Galactocentric circles with radii of 2 to 18 kpc 
are overplotted on the map. 

By adopting the entire sample we calculated once again the iron Galactic gradient. Note that 
the new sample includes 60 Cepheids located inside the Solar circle together with two Cepheids 
located beyond 16 kpc. We found that the slope of the linear iron gradient ov er the w hole sam- 
ple is -0.052+0.003 dex kpc ' (Fig. |5]l. This slope is slightly s hallower than in IPaper j where we 
found -0.061+0.019 dex k pc~', as we l l as in 



Friel et al 



Luck et al 



(l2002h who found -0.06+0.01 dex kpc" 



from open clusters and in 
Cepheids. Data plotted in this figure indicate that: 



(l2006l) who found -0.068+0.003 dex kpc ' from classical 



- The spread in metalhcity is larger in the outer disk, in particular around 10 kpc. 

- There is no evidence of a gap in the Galactocentric distances covered by the current sample. 
We will discuss this point in the next section. 
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5 10 -5 20 

X-AxIs of Galactocentric Cartesian Coordinates (kpc) 

Fig. 4. Spatial distribution of the entire Cepheid sample. Filled diamonds mark our data, while open 
circles show data from Andrievsky et al. and open triangles are data from Mottini et al. The x-axis 
contains both the center of the Galaxy at (0,0) and the Sun at (8.5,0). Galactocentric circles with 
radii of 2 to 18 kpc are overplotted on the map. 



1 .0 



0.5 



_ SLOPE: -0.052 dex/kpc 



0.0 



-0.5 



-1.0 




8 10 12 14 

Galactocentric distance ( kpc ) 



1 6 



Fig. 5. Galactic iron abundance gradient. Filled diamonds mark our data, open circles show data 
from Andrievsky et al. and open triangles data from Mottini et al. The solid line shows the linear 
regression over the entire sample. Characteristic error bars are plotted in the lower right corner. 



4.2. The Iron gradient in the outer disk 

To properly investigate the iron gradient in the outer disk, we selected (in our sample of 63 stars) 
Cepheids located between 10 and 15 kpc (see Fig.|6]). We estimated the Galactic gradient and found 
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Table 5. Additional stars from Andrievsky's series and from 
are marked with an (*). 



Romaniello et al 



(l2008h . The latter 



Distances are based on infrared photometry and [Fe/H] found in the literature. 



Object 




[Fe/H] 


Object 




[Fe/H] 


Object 


d;ft 


[Fe/H] 


Object 




[Fe/H] 


Object 




[Fe/H] 




kpc 


dex 




kpc 


dex 




kpc 


dex 




kpc 


dex 




kpc 


dex 


TT Aql 


7.70 


0.11 


BD Cas 


9.85 


-0.07 


MW Cyg 


8.15 


0.09 


S Nor 


7.74 


0.05 


RY Sco 


7.16 


0.09 


FM Aql 


7.86 


0.08 


CFCas 


10.30 


-0.01 


V386 Cyg 


8.48 


0.11 


V340 Nor 


6.90 


0.00 


KQ Sco 


6.26 


0.16 


FN Aql 


7.41 


-0.02 


CH Cas 


9.98 


0.17 


V402 Cyg 


8.18 


0.02 


Y Oph 


7.94 


0.05 


V500 Sco 


7.08 


-0.02 


V496 Aql 


7.63 


0.05 


CY Cas 


10.48 


0.06 


V532 Cyg 


8.55 


-0.04 


BE Oph 


7.69 


0.00 


SS Set 


7.61 


0.06 


VI 162 Aql 


7.48 


0.01 


DD Cas 


10.22 


0.10 


V924 Cyg 


7.92 


-0.09 


SV Per 


10.82 


0.01 


UZ Set 


5.69 


0.33 


Eta Aql 


8.30 


0.05 


DF Cas 


10.45 


0.13 


VI 334 Cyg 


8.46 


-0.04 


UXPer 


12.30 


-0.21 


EW Set 


8.14 


0.04 


V340 Ara 


5.00 


0.31 


DL Cas 


9.46 


-0.01 


VI 726 Cyg 


8.73 


-0.02 


VX Per 


10.46 


-0.05 


V367 Set 


6.81 


-0.01 


RT Aui- 


8.93 


0.06 


FM Cas 


9.58 


-0.09 


Beta Dor 


8.50 


-0.01 


AS Per 


9.79 


0.10 


BQ Ser 


7.70 


-0.04 


RX Aui- 


10.03 


-0.07 


V636 Cas 


8.96 


0.06 


WGem 


9.44 


-0.04 


AW Per 


9.29 


0.01 


SZ Tau 


8.94 


0.08 


AN Aur 


11.69 


-0.16 


V Cen 


7.99 


0.04 


Zeta Gem 


8.86 


0.04 


BM Per 


11.30 


0.10 


AE Tau 


12.33 


-0.19 


CY Aui- 


13.58 


-0.40 


CP Cep 


9.79 


-0.01 


V Lac 


9.17 


0.00 


HQ Per 


13.95 


-0.31 


S Vul 


7.60 


-0.02 


ER Aui- 


16.69 


-0.34 


CR Cep 


9.00 


-0.06 


XLac 


9.01 


-0.02 


V440 Per 


9.13 


-0.05 


T Vul 


8.35 


0.01 


RW Cam 


10.01 


0.04 


IR Cep 


8.65 


0.11 


Y Lac 


9.08 


-0.09 


X Pup 


10.50 


-0.03 


U Vul 


8.16 


0.05 


RXCam 


9.25 


0.03 


Del Cep 


8.57 


0.06 


Z Lac 


9.20 


0.01 


S Sge 


8.13 


0.10 


X Vul 


8.10 


0.08 


TV Cam 


12.38 


-0.08 


BG Cru 


8.30 


-0.02 


RRLac 


9.21 


0.13 


USgr 


7.89 


0.04 


SV Vul 


7.80 


0.03 


AB Cam 


12.72 


-0.09 


XCyg 


8.32 


0.12 


BG Lac 


8.77 


-0.01 


WSgr 


8.08 


-0.01 


U Car 


8.11 


0.12* 


AD Cam 


13.47 


-0.22 


SU Cyg 


8.17 


0.00 


V473 Lyr 


8.33 


-0.06 


YSgr 


8.02 


0.06 


WZCar 


8.14 


0.18* 


RWCas 


10.39 


0.22 


SZ Cyg 


8.60 


0.09 


TMon 


9.75 


0.22 


VY Sgr 


6.29 


0.26 


VW Cen 


6.83 


0.05* 


RYCas 


9.91 


0.26 


TXCyg 


8.46 


0.20 


AA Mon 


12.03 


-0.21 


WZ Sgr 


6.71 


0.17 


XX Cen 


7.52 


0.04* 


SU Cas 


8.71 


-0.01 


VXCyg 


8.58 


0.09 


CU Mon 


14.73 


-0.26 


XX Sgr 


7.18 


0.10 


KNCen 


6.92 


0.07* 


SW Cas 


9.34 


0.13 


VY Cyg 


8.49 


0.00 


EE Mon 


16.06 


-0.51 


APSgr 


7.68 


0.10 


GH Lup 


7.56 


0.03* 


SY Cas 


9.66 


0.04 


VZCyg 


8.75 


0.05 


EG Mon 


14.43 


-0.20 


AV Sgr 


6.37 


0.34 


S Mus 


8.12 


0.18* 


SZ Cas 


10.25 


0.04 


BZ Cyg 


8.54 


0.19 


EI Mon 


13.12 


-0.18 


BB Sgr 


7.71 


0.08 


UU Mus 


7.60 


0.05* 


TU Cas 


8.95 


0.03 


CD Cyg 


8.05 


0.07 


V504 Mon 


11.59 


-0.31 


V350 Sgr 


7.62 


0.18 


UNor 


7.34 


0.10* 


XY Cas 


9.73 


0.03 


DTCyg 


8.42 


0.11 


V526 Mon 


9.69 


-0.13 


RV Sco 


7.74 


0.10 


LS Pup 


11.39 


-0.11' 



a slope of -0.012+0.014 dex kpc"'. The shallow slope clearly indicates a flattening of the iron gra- 
dient in the outer disk. This also suggests that the Galactic iron gradient might be more accurately 
described by a bimodal distri bution. This evidenc e agrees quite well with the slope (-0.018+0.02 
dex kpc"') recently found bv lCarraro et alj (l2007h using old Open Clusters ranging from 12 to 21 
kpc and a Galactocentric distance of the Sun of 8 .5 kpc. It is also similar to the flat distribution 
(+0.004+0.01 1) found bv Undrievskv et al 



(120041) in the outer zone (zone III) of their multimodal 
gradient. Moreover and even more importantly, the two quoted independent investigations suggest, 
within the errors, a similar result concerning the value of the metallicity ba sement. The value of 
the metallicity basement we found from our data is -0. 19+0. 14 dex, whereas ICarraro et al.l (l2007b 
found that open clusters located in the outer disk have an iron abu ndance of [Fe/H] - -0.34+0.15 
dex, a value close to the — 0.3 dex found bv lTwarog et al. ( 1997) from the sam e kind of objects. 
Our basement value also agrees very well with the one of lAndrievskv et al.l (120041) as the mean iron 
abundance obtained in their outer disk zone is -0.209+0.058 dex. 

Current findings a re marginally consistent with the results concerning the metallicity basement 
bv lYong et al.l (l2006h . They suggested a two-level basement in the outer disk: one level located 

at 0.5 dex, is notably lower than the other values, while the second one is still lower at 

0.8 dex. According to these authors, the more metal-poor values might be explained by Cepheids 
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Fig. 6. Same as Fig.|3] but for Cepheids with Galactocentric distances ranging from 10 to 15 kpc. 



acquired by the Galaxy through merging events. This working hypothesis was also suppo r ted by 
the significantly differen t [a/Fe] ratios shown by these objects. In this context. I Yong et al.l ( l2005h 
and 



Carney et al 



(l2005h using different stellar tracers, namely open clusters and field red giants, 
found a metallicity basement of ^ -0.5 dex. 

If we now study the same region including all the Cepheids for which we have accurate esti- 
mates of both iron abundance and distance, the gradient is affected by selection effects. The value 
of the slope depends on: 

- the range of Galactocentric distances covered by the Cepheid sample, since the targets are not 
uniformly distributed along the disk. 

- the spatial sampling across the disk. By using Cepheids in our sample with R > 10 kpc we 
found a slope of- -0.012±0.014 dex kpc On the other hand, if we use Cepheids from the 
Andrievsky sample and our distances we found a slope of -0.077+0.010 dex kpc while with 
two samples together we obtain a slope of -0.050+0.008 dex kpc"' (See Fig. Ell. 

The quoted results are strongly affected by the poor spatial sampling. The number of known 
Cepheids located in the outer disk is quite Umited and quite often they are also relatively faint. This 
means that long exposure times with 4-8m class telescopes are necessary to obtain the S/N required 
for accurate spectroscopy. 

The other interesting findings brought forward by the data plotted in Fig.|2]is the large spread in 
metallicity between 10 and 12 kpc. This feature was partially expected and is caused by the decrease 
in gas density and in star formation rate when moving from the innermost to the outermost regions 
of the Galaxy. 

- A group of Cep heids has rather hig h metallicities ( — hO.O to — 1-0.2 dex). Some of them were 
already noted by|] 



Luck et al. 



(12006 1) to be located approximately at the same Galactic longitude 



120' 
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Fig. 7. Galactic iron abundance gradient. Filled diamonds mark our data, while open circles show 
data from Andrievsky et al. and open triangles display data from Mottini et al. The solid line 
shows the slope based on our sample, the dotted line the slope based on the additional data 
from Andrievsky et al. and the dashed line the slope obtained using the two quoted samples. 
Characteristic error bars are plotted in the lower right corner. 



- On the contrary, a group of stars has rather low metallicities ( 0.2 to 

them lying in the 180°< I < 200°range of Galactic longitudes. 



-0.4 dex), most of 



According to this circumstantial empirical evidence, if we plot the Cepheids located between 
10 and 12 kpc on a Galactocentric polar map (see Fig.[8]l, we can see that the metallicity is de- 
pendent on the Galactocentric longitude, with a spread in metallicity reaching 0.8 dex for stars 
located at roughly the same Galactocentric distance. Here we define the Galactocentric longitude 
as the angle of a polar coordinate system centered on the Galactic center. Galactocentric longitudes 
are measured from -180°to +180°, counter clockwise from the origin, arbitrarily defined as the 
Galactic radius containing the Galactic center at (0,0) and the Sun at (8.5,0). Distances projected 
on Galactocentric Cartesian coordinates are also labelled on the x and y-axis. 

The occurrence of a large spread in metal abundance between stars located at the same radius 
of the Galactic disk allows us to explain why we do not observe a gap in the Galactic gradient 
as proposed in some previous studies. The gap might be caused either by selection effects in the 
sampling of the disk, or by local inhomogeneities in the disk, and therefore, restricted to a narrow 
range in longitude. The C epheid sample collected by 



Andrievsky et al. 



Andrievskv et al 



( I2002cl) : 



Luck et al. 



(I2003h : 



(120041) are indeed located in a limited range of Galactic longitu de (190°< / < 



250°). The same outcome applies to the sample of 76 Open Clusters collected by iTwarog et al 
119971) (130°</< 260°). 



Current findings suggest the need to be very cautious in using the slopes of Galactic gradients, 
in particular in the outer disk due to poor spatial sampling. A significant improvement in the sample 
of outer disk Cepheids requires a systematic photometric survey to detect new Cepheids. This goal 
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Fig. 8. Metallicity of Cepheids located between 10 and 12 kpc. Filled diamonds mark our data 
while filled circles show data from Andrievsky et al. The metallicity strongly depends on the 
Galactocentric longitude (see text). The center of the Galaxy is at (0,0) and the Sun is at (8.5,0). 
Galactocentric circles with radii of respectively at 6, 8, 10, 12 and 14 kpc are over-plotted on the 
map. 



can be achieved with the new generation of automatic telescopes such as HAT (Hungarian-made 
Automated Telescope; http ://cf a- w w w.harvard . edu/~gbakos/H AT/index .html) ! or ASAS (All Sky 
Automated Survey; http://archive.princeton.edu/~asas/l l and with spectroscopic follow up of the 
new targets. 



5. Conclusions 

High resolution spectra and accurate distance determinations based on infrared photometry allowed 
us to extend our study of the Galactic iron abundance gradient toward the outer disk. New data 
led us to a slope of -0.052+0.003 dex kpc ' in the 5-17 kpc range, in very good agreement with 
previous studies if we consider a linear gradient. However, the Galactic iron abundance gradient 
is better described by a bimodal distribution with a steeper slope toward the bulge and a flattening 
of the gradient toward the outer disk. In particular, we found a shallow slope of -0.012+0.014 
dex kpc ' in this region (10-15 kpc). CuiTent data show no evidence of a sharp discontinuity in 
metallicity for Galactocentric radii of ~ 10-12 kpc. Poor spatial sampling aff'ecting the adopted 
stellar tracers do not allow us to constrain on a quantitative basis whether the occurrence of this 
gap is real or caused by local inhomogeneities in the metallicity distribution of the disk. We also 
found a large spread in metal abundance (« 0.8 dex) for Cepheids located approximately at the 
same distance (10 < R < 12 kpc), but covering a broad range of Galactocentric longitude. A larger 
Cepheid sample, uniformly distributed over the four Galactic quadrants, is mandatory to constrain 
this eff'ect. 
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